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£j<act  AlgorithMs 

In  this  section  it  will  be?  assumed  that  the  letter  n  is 
a  positive  integer.  If  l<  e  C.0 , 1 ,2,  .  .  .  ,n3 ,  let  B^  n  denote  the 
graph  on  n  vertices  whose  vertices  are  numbered  Ct,2,...,n> 
and  whose  edge  set  is  CCa  ,b}  ;  la  -b  Kk5  .  Given  a  graph  G  on  n 
vertices  we  can  determine  its  bandwidth  as  follows 
(See  CIO! . ) : 

Step  1.  Set  m  <=  0 

Step  2.  If  G  is  a  subgraph  of  Bra  n  then  stop,  m 
is  the  bandwidth 

Step  3.  Set  m  (~  m  +  1  and  return  to  Step  2 


Unfortunately,  this  algorithm  has  little  practical  sig¬ 
nificance  since  determining  whether  or  not  a  graph  is  a  sub¬ 
graph  of  another  graph  takes  a  prohibitive  amount  of  time. 

In  1976,  Papad imi t r i ou  CIO]  was  able  to  prove  that  the  band¬ 
width  minimization  problem  is  (in  general)  NP-complete.  In 
1978,  Garey,  Graham,  Johnson,  and  Knuth  C7]  were  able  to 
prove  NP-completeness  in  a  sharper  form.  In  the  same  paper 
they  give  a  linear  algorithm  for  determining  whether  or  not  a 
graph  has  bandwidth  2.  Gibbs  and  Poole  CIO)  and  Fulkerson 
and  Gross  16]  had  previously  presented  linear  time  algorithms 
for  bandwidth  1.  The  only  gap  now  is  the  question  of  whether 
or  not  the  problem  of; 
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Is  bandwidth(G)  <  3  NP-complete  for  arbitrary  G? 


Pr  od  uc  t  i  on  Algor  ithfjs 

Large  sparse  Matrix  problems  arise  in  a  variety  of 

application  areas,  such  as  structural  engineering,  fluid 

dynamics,  and  network  analysis.  Although  out  of  vogue  for  a 

while,  band  solvers  which  solve  sparse  linear  systems  of 

equations  directly  using  Gaussian  elimination  have  been  used 

increasingly  on  vector  processing  Machines  such  as  the  Cray 

and  the  CDC  Star.  <See  [153.)  Pure  Gaussian  elimination 
3 

requires  o(n  )  operations  for  an  n  x  n  matrix.  The  time 

2 

required  to  solve  the  systeM  using  a  band  solver  is  o(nb  ) 
where  b  is  the  bandwidth.  <If  A  is  an  n  x  n  Matrix  its  band¬ 
width  is  Max  li-il.)  The  probleM  then  becones  how  to  reduce 

V° 

b  as  much  as  possible.  Apparently  people  bee a Me  adept  at 
labeling  by  hand  the  grids  (graphs)  which  the  Matrices  repre¬ 
sented.  The  motivation  then  existed  for  the  creation  of 
automatic  production  programs  which  (hopefully)  reduced  band¬ 
width. 


In  a  1979  survey  paper  by  Everstine  [53,  49  "reduction" 
algorithms  are  referenced.  This  list  is  by  no  means  exhaus¬ 
tive  and  it  is  safe  to  assume  that  between  50  and  100  reduc¬ 
tion  algorithms  h.v»e  been  designed  to  date.  Although  most 
applications  pronrammers  working  in  the  area  were  not  famil- 
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i ar  with  the  results  on  NP-completeness,  hueristic  reduction 
algorithms  were  developed.  Most  were  designed  with  the  goal 
of  obtaining  acceptable  (rather  than  exact)  bandwidth. 

Alway  and  Martin  published  the  first  bandwidth  reduction 
algorithm  in  1965  [13.  The  Main  idea  of  the  algorithm  was  to 
examine  Many  row  and  corresponding  coluwn  perMutations  with 
sone  rather  conplex  criteria  used  to  discard  certain  perMuta¬ 
tions.  The  authors  adMit  in  the  paper  that  the  algorithm 
takes  too  Much  tiMe  even  for  graphs  of  Moderate  size. 
Nevertheless,  they  created  an  algorithm  which  was  effective 
for  their  use  in  reducing  bandwidth  of  snail  Matrices  and 
stiMulated  further  research. 

The  first  widely  used  production  algorithm  was  published 
by  Rosen  [193  in  1968.  The  idea  behind  the  algorithm  is  to 
first  coMpute  the  bandwidth  of  the  Matrix  (graph).  The  end¬ 
points  of  the  edges  which  cause  the  bandwidth  to  be  attained 
are  then  exaMined  to  see  if  their  vertex  labels  can  be  inter¬ 
changed  with  other  vertex  labels.  When  such  row  and  corres¬ 
ponding  coluMn  pernutations  are  performed,  the  bandwidth  is 
reconputed  to  see  if  an  improvement  was  Made.  The  entire 
process  iterates  until  no  improvement  to  bandwidth  is  Made. 
The  algorithm  uses  a  local  strategy  in  terms  of  the  entire 
structure  of  the  associated  graph,  and  is  tiMe  consuming. 

This  algorithm  was  published  along  with  its  imp lemen t a t i on  as 
a  FORTRAN  progr.m.  Perhaps  the  Most  interesting  aspect  of 
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the  algorithm  is  that  people  were  still  using  it  in 
production  Mode  in  the  late  1970's. 

A  Major  breakthrough  occurred  in  1969  when  the  faMous 
Cuthi 11-McKee  algorithm  [31  was  published  and  subsequently 
incorporated  into  hiASTRAN  <See  [43.).  To  understand  the  Cut- 
hill-McKee  algorithm  and  the  algorithm  which  follows,  it  is 
necessary  to  define  a  few  terMS.  The  definitions  and  the 
next  two  algorithm  descriptions  are  taken  fron  [111.  A  level 
structure,  L<G>,  of  a  graph  G  is  a  partition  of  the  vertices 
of  the  graph  into  levels  ,L2  ,...,Lksuch  that 

1.  all  vertices  adjacent  to  vertices  in  level  , 
are  in  either  level  L^or  L2, 

2.  all  vertices  adjacent  to  vertices  in  level  l_k 
are  in  either  level  L^or  Lj^ 

and 

3.  for  l<i<k,  all  vertices  adjacent  to  vertices 
in  level  L^are  in  either  level  Li_1  ,1^,  or  L_i+1  . 

To  each)  vertex  v  of  the  graph  there  corresponds  a  parti¬ 
cular  level  structure  L  < G >  called  the  level  structure  rooted 
at  v.  Its  levels  are  determined  by 

1  .  1-2=  ( vl  ,  and 

2.  for  i>l,  L^is  the  set  of  all  those  vertices  which  are 
not  yet  assigned  to  a  level,  hut  are  adjacent  to  ver- 
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tices  of  level  L  ^ 

In  any  level  structure  L(G)  ,  rooted  or  not,  w^  <L>  = 

I  L.  I  is  called  the  width  of  level  i,  and  w(L)  =  MaxCw.l  is 
~i  -  ~~  . . .  i  1 

the  width  of  the  level  structure  L(G).  It  is  easily  observed 

that  for  any  level  structure,  L,  a  numbering  of  G  which 

assigns  consecutive  integers  level  by  level,  first  to  the 

vertices  of  level  L  ,  then  to  those  of  level  L  ,  and  so 

“1  *  ~2 

forth,  gives  a  bandwidth  which  is  less  than  or  equal  to 
2w(L)-l ,  Furthermore,  if  the  level  structure  is  rooted,  then 
the  bandwidth  of  the  numbering  is  greater  than  or  equal  to 
w(L) .  The  depth  of  a  level  structure  is  the  number  of  lev¬ 
els,  k. 

For  discussion  purposes  assume  all  graphs  are  connected, 
The  Cuthill-McKee  algorithm  begins  by  generating  a  level 
structure  rooted  at  each  vertex  of  low  degree.  See  131  for 
their  definition  of  low  degree.  For  each  rooted  level  struc¬ 
ture  of  Minimal  width  generated  in  the  first  step  the  graph 
is  numbered  level  by  level,  with  consecutive  positive  integers 
by  first  assigning  the  integer  1  to  the  root,  For  each  suc¬ 
cessive  level,  starting  with  level  2,  the  vertices  adjacent 
to  the  lowest  numbered  vertex  in  the  previous  level  are  num¬ 
bered  according  to  increasing  degree  with  ties  broken  arbi¬ 
trarily.  The  remaining  vertices  adjacent  to  the  next  lower 
numbered  vertex  of  the  preceding  level  are  numbered  next, 
again  in  order  of  increasing  degree.  The  process  continues 
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until  all  vertices  of  the  current  level  are  numbered.  This 
numbering  process  is  repeatedly  applied  in  order  to  subse¬ 
quent  levels  until  all  vertices  have  been  numbered.  The  num¬ 
bering  which  gives  the  best  bandwidth  of  any  of  the  rooted 
level  structures  of  Minimal  width  is  selected.  This  algor¬ 
ithm  was  the  most  widely  used  bandwidth  reduction  algorithm 
during  the  1970's. 


In  1976  Gibbs,  Poole,  and  Stockmeyer  published  a  series 
of  related  papers  <See  till,  [21,  and  [121.)  describing  a  new 
bandwidth  reduction  algorithm  which  overcame  what  they  per¬ 
ceived  to  be  major  shortcomings  of  the  Cuthill-McKee  algor¬ 
ithm.  The  algorithm  has  three  phases,  finding  a  pseudo-diam¬ 
eter,  minimizing  level  width,  and  numbering  a  level  structure 
(not  necessarily  rooted). 

The  Cuthill-McKee  algorithm  is  inefficient  because  of 
the  time  consumed  performing  an  exhaustive  search  to  find 
rooted  level  structures  of  minimal  width.  A  second  problem 
is  that  the  graph  is  renumbered  and  the  corresponding  band¬ 
width  recomputed,  for  every  level  structure  found  of  Minimal 
width.  A  third  problem  is  that  the  bandwidth  obtained  by  a 
Cuthill-McKee  numbering  can  never  be  less  than  the  width  of 
the  rooted  level  structure  used,  although  the  bandwidth  of  a 
graph  can  be  less  than  the  width  of  any  rooted  level  struc¬ 
ture.  The  first  two  shortcomings  are  overcome  by  carefully 
selecting  a  starting  vertex  after  generating  only  a  rela- 
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tively  SMall  n uMher  of  level  structures,  The  graph  is 
renutibered,  and  corresponding  bandwidth  conputed  only  once. 
The  third  problem  is  resolved  by  utilizing  a  wore  general 
type  of  level  structure, 

Gibbs,  Poole,  and  StockMeyer  observed  that  level  struc¬ 
tures  of  sMall  width  are  usually  among  those  of  maximal 
depth.  Clearly,  increasing  the  number  of  levels  always 
decreases  the  average  number  of  vertices  in  each  level,  and 
tends  to  reduce  the  width  of  the  level  structure  as  well. 
Ideally,  then  one  would  like  to  generate  level  structures 
rooted  at  endpoints  of  a  diameter.  Since  there  is  no  known 
efficient  procedure  which  always  finds  such  vertices,  the 
following  algorithm  is  used  to  find  the  endpoints  of  a  R§eu- 
dfizdiameter ,  that  is,  a  pair  of  vertices  that  are  at  nearly 
maximal  distance  apart.  For  a  large  class  of  graphs,  includ¬ 
ing  all  trees,  the  pseudo-diameter  is  actually  a  real  diame¬ 
ter  . 


Finding  the  endpoints  of  a  pseudo-diaMeter 

A.  Pick  the  smallest  numbered  vertex  of  minimal  degree  and 
call  it  v . 

B.  Generate  a  level  structure  l_v  rooted  at  vertex  v.  Let  S 
be  the  set  of  vertices  which  are  in  the  last  level  of  L 
(i.e.,  those  vertices  which  are  furthest  away  from  v). 

C.  Generate  level  structures  rooted  at  vertices  se  S 
selected  in  order  of  increasing  degree.  If  for  some  seS  the 
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depth  of  L_s  is  greater  than  the  depth  of  l_v  ,  then  set  v  <  =s 
and  return  to  step  B. 

D,  Let  u  be  the  vertex  of  S  whose  associated  level  structure 
has  shii  1  lest  width.  The  algorithM  terwinates  with  u  and  v 
the  endpoints  of  a  pseudo- diafteter . 

In  the  process  of  finding  a  pseud  o-d:i.aMeter ,  this  algor¬ 
ithM  constructs  level  structures  L  and  L  rooted  at  the  end- 

u  v 

points  u  and  v,  respectively.  It  is  possible  to  cowbine 
these  two  level  structures  into  a  new  level  structure  whose 
width  is  usually  less  than  that  of  either  of  the  original 
rooted  ones,  using  the  following  algorithM. 

MiniMizing  Level  Width 

A.  Using  the  rooted  level  structures  Ly  =  CL^  >l_2 ,  .  .  .  ,Lk3  and 
Lu  -  ,M2 ,  .  .  .  obtained  froM  the  algo.rithM  for  finding  a 

pseud o-d iaMeter ,  associate  with  each  vertex  w  of  G  the  ord¬ 
ered  pair  <i,j),  called  the  associated  level  pair,  where  i  is 
the  index  of  the  level  in  Ly  that  contains  w,  and  k+l-j  is 

the  index  of  the  level  in  L  that  contains  w.  Thus  the  pair 

u 

<i,j)  is  associated  with  a  vertex  w  if  and  only  if 
w  e  L.  j  ‘  Note  that  the  pair  (1,1)  is  associated  with 

the  vertex  v,  while  the  pair  ( k ,  l<  >  is  associated  with  u. 

E< .  Assign  the  vertices  of  G  to  levels  in  a  new  level  struc¬ 
ture  L.  ~  ,N2  ,  .  .  .  ,N^3  as  follows: 

1.  If  the  associated  level  pair  of  a  vertex  w  is  of  the 
forM  <i,i>,  then  vertex  w  is  placed  in  N,  . 


The  vertex 
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w  and  all  edges  incident  to  w  are  renoved  Iron  the 
graph.  If  V(G)  is  empty,  stop. 

2.  The  graph  G  now  consists  of  a  set  of  one  or  More 
disjoint  connected  components  C^,C2>.".>Ct  ordered  so 

that  i v < >  mv<c2>  \  >_.  .  <2,v(Ct  >  1  • 

3.  For  each  connected  component  (  taken  in  the  order  , 
C2>...>Ct)  do  the  following: 

a.  Compute  the  vector  <  n-^  ,  >  •  •  •  >  >  where 

nt  »  I  Hi  I  . 

b.  Compute  the  vectors  < h^ , h2 >  •  •  •  > h^ >  and 

<  1^  ,  12  >  .  .  .  >  1^  )  where  -  n  +  (the  number  of  ver¬ 
tices  which  would  be  placed  in  N.  if  the  first 
element  of  the  associated  level  pairs  were  used) 
and  li  =  ni  +  (the  number  of  vertices  which  would 
be  placed  in  ^  .if  the  second  element  of  the  asso¬ 
ciated  level  pairs  were  used). 

c.  Find  hg  =  max  ih^-n^  >03  and 
ln  -  max  C 1 . : 1 . -n . > 03 . 

U  111 

i.  If  h0<l0,  place  all  the  vertices  of  the 
connected  component  in  the  levels  indicated  by 
the  first  elements  of  the  associated  level 
pairs. 

ii.  If  10  <  h0  ,  use  the  second  elements  of  the 
level  pairs  to  place  the  vertices  in  the 
levels. 

iii.  If  hg  ~1q  then  use  the  elements  of  the 

level  pairs  which  arise  from  the  rooted  level 
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structure  of  sMallor  width.  If  the  widths  are 
equal,  use  the  first  elements. 

The  algorithM  ter ni nates  when  each  vet ex  of  G 
has  been  assigned  a  level  in  the  level  structure  L. 

The  numbering  procedure  is  similar  to  that  of  the  Cut- 
hill-McKee  algorithM  in  that  is  assigns  consecutive  positive 
integers  to  the  vertices  of  G  level  by  level.  A  few  Modifi¬ 
cations  were  necessary,  however,  since  the  level  structures 
obtained  by  the  algorithM  to  Minimize  level  width  are  of  a 
More  general  type  than  the  rooted  ones  used  by  the  Cuthill- 
McKee  algorithM.  Under  certain  conditions,  profile  can  be 
further  reduced  by  using  the  reverse  nunbering  described  in 
step  D  below.  (See  the  next  section  on  profile  reduction.) 

NuMber ing 

A.  If  the  degree  of  u  is  less  than  the  degree  of  v,  then 

interchange  o  and  v  and  reverse  the  level  structure  obtained 

by  the  algorihM  to  Minimize  level  width  by  setting  N  to 

i 

N  .  (This  insures  that  the  nuMhering  starts  froM  the  end- 

k-i+1 

poirit  of  lower  degree). 

B.  Assign  consecutive  positive  integers  to  the  vertices  of 
level  in  the  following  order: 

1.  Assign  the  nuMber  1  to  the  vertex  v. 

2.  Let  w  be  the  lowest  nur.bered  vertex  of  level 
which  has  unnumbered  vertices  in  N  adjacent  to  it. 
Nunher  the  vertices  of  adjacent  to  w,  in  order  of 
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incr e a s i n q  d e q rcc  ,  R e p e a t  t h i s  s t ep  until  all  vertices 
of  N  adjacent  to  numbered  vertices  are  thenselves  nuM- 
bered . 

3.  If  any  unnunbered  vertices  renain  in  level  ,  nuh- 
ber  the  one  of  Minimal  degree,  then  go  to  step  B.2. 
Otherwise  proceed  to  step  C. 

C.  Number  the  vertices  of  level  N  ,  i~2,3, . . . ,k  as  follows: 

1.  Let  w  be  the  lowest  numbered  vertex  of  level  N.  , 

i-I 

that  has  unnumbered  vertices  of  level  N  adjacent  to  it. 

Nunber  the  vertices  of  N  adjacent  to  w  in  order  of 

increasing  degree.  Repeat  this  step  until  all  vertices 

of  level  N.  adjacent  to  vertices  of  level  N.  ,  are  num- 
i  i-I 

bered . 

2 .  R ep  ea t  s t  ep s  B .2  and  B .  3 ,  r  ep 1 a  c i n  g  1  with  i . 

D.  The  nunberinq  is  reversed  by  setting  :i  to  n-i  +  1  ,  for 

i--l  , 2 >  .  .  .  ,n  if  either  of  the  two  following  conditions  holds: 

1.  Step  A  interchanged  vertices  u  and  v  and  the  algor- 
ithn  for  minimizing  level  width  selected  the  second  ele~ 
Ments  of  the  level  pairs  for  component  C  . 

2.  Step  A  did  not  interchange  vertices  u  and  v  and  the 
algorithm  for  Min  ini  zing  level  width  ..elected  the  first 
elements  of  the  level  pairs  for  component  . 

Recently  Lewis  f  151  has  made  some  improvements  to  the 
FORTRAN  implementation  of  the  Gibbs,  Poole,  and  Stockmeyer 
algorithm.  He  claims  that  this  new  version  will  run  slightly 
faster  than  Algorithm  50C  [21. 
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Prof il e  R e d u c t i o n 

A  related  but  independent  problem  is  that  of  profile 

reduction.  Let  A  be  an  n  x  n  sparse  symmetric  Matrix.  To 

define  the  profile  of  A,  first  define  f  =  ninCjia  •/  03  for 

i  -* 

i~l,2,...,n  (assuMe  aii^  0).  This  locates  the  leftmost 
nonzero  element  in  each  row.  Next  define  p.=  i-f^.  The  pro¬ 
file  is  defined  to  be 

n 

E  Pi 
i-1  1  * 

Note  that  in  the  figure  below  the  first  labeling  gives  a 
bandwidth  of  3  (which  is  Minimal)  and  a  profile  of  8.  The 
second  labeling  gives  a  bandwidth  of  5  and  a  profile  of  5 
(which  is  Minimal). 


The  first  practical  and  widely  used  profile  reduction 
algorithm  was  the  reverse  Cuthill-McKee  algorithm  which  was  a 
modification  of  the  original  algorithm  due  to  George  181. 

The  algorithm  is  basically  the  same  as  described  above  except 
a  step  is  added  at  the  end  which  "reverses"  the  numbering  by 
replacing  i  with  n-i+1 .  The  reversal  of  the  numbering  has  no 
effect  on  bandwidth  and  Liu  and  Sherman  Cl 71  have  proved  that 
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it  can  newer  increase  profile. 

Other  good  profile  reduction  algorithms  have  been  devel¬ 
oped  by  King  114.1,  Gibbs  [133,  and  On  ay  [20  3.  Snay's  paper 
contains  an  exellent  description  of  the  intuitive  notions  of 
why  the  King  and  Cu t h i 1 1 -McKee  numberings  are  good  for  pro¬ 
file  and  bandwidth  reduction  r espec t i vel y .  Of  the  profile 
reduction  algorithms  Mentioned,  Snay's  tends  to  produce  the 
best  (least)  profile  when  it  works  well.  As  will  be  Men¬ 
tioned  below,  Snay's  algorithm  sometimes  has  trouble  as  does 
King's  on  certain  graphs.  Recently,  Lewis  [153  modified  the 
FORTRAN  implementation  of  Gibbs'  profile  reduction  algorithm 
so  that  it  runs  almost  as  fast  as  the  Gibbs,  Poole,  and 
Stockmeyer  bandwidth  reduction  algorithm.  Although  Gibbs' 
profile  reduction  algorithm  does  not  always  reduce  profile  as 
much  as  Snay's  algor ithm, Lewis'  implementation  is  a  signifi¬ 
cant  contribution  in  that  it  has  Made  the  More  stable  algor¬ 
ithm  fast  enough  to  be  a  practical  alternative  to  Snay's 
algorithm  and  King's  algorithm.  Based  on  percentage  of 
reduction,  one  can  claim  that  the  profile  reduction  of  the 
three  algorithms  are  comparable. 


Evaluation  of  Reduction  Algorithms 


Since  bandwidth  and  profile  reduction  algorithms  are 
hueristic,  it  is  difficult  to  compare,  them.  An  approach 
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taken  by  Gibbs,  Poole,  and  StockMeyer  112.1  has  gained  in 
popularity  and  has  probably  defined  the  de  facto  standard  for 
reduction  algorithm  testing.  Dr.  Gordon  Everstine  of  the 
David  Taylor  Naval  Ship  and  Research  Development  Center  in 
Bethesda,  Maryland  (DTNSRDC)  provided  then  with  19  test 
Matrices  which  arise  when  the  finite  element  Method  is  used 
to  approximate  solutions  of  partial  differential  equations  in 
several  structural  engineering  applications.  The  structures 
included  parts  of  aircraft,  gasoline  storage  tanks,  subma- 
rines,  propellor  blades,  and  satellites.  More  recently, 
Everstine  [51  asseMbled  a  new  set  of  application  Matrices. 
This  set  of  30  Matrices  includes  some  of  the  original  19  and 
drawings  of  the  structutes  represented  by  the  Matrices  are 
available.  This  set  of  Matrices  is  now  the  standard  bench- 
Mark  for  any  production  reduction  algorithm. 

In  the  sane  paper,  Gibbs,  Poole,  and  Stockmeyer  ran  a 
second  set  of  tests.  They  were  interested  in  obtaining  tim- 
ing  data  in  order  to  estimate  the  asyMptotic  behavior  of  the 
algorithms  they  were  testing.  The  algorithms  were  tested  on 
several  families  of  grids.  The  grids  selected  were  n  x  n 
squares,  3r>  x  n  rectangles,  n  x  20n  rectangles,  20  rt  x  n  cyl¬ 
inders,  and  n  x  20n  cylinders.  (Bee  the  paper  for  More 
details  about  the  types  of  elements  used  in  the  families  of 
grids.)  George  and  Liu  adopted  the  saMe  technique  in  Liu's 
thesis  t 1 63  and  in  a  later  paper  193. 
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The  problem  of  how  to  evaluate  a  heuristic  algorithm  has 
no  easy  solution,  It  is  usually  not  difficult  to  cone  up 
with  bizarre  c o un ter examp ies  which  will  destroy  the  perfor¬ 
mance  of  almost  any  reduction  algorithm  ~~  good  or  bad. 
Experience  has  shown,  however,  that  the  results  of  testing 
algorithms  on  the  set  of  19  test  matrices  and  then  on  the  set 
of  30  test  matrices  correlate  well,  These  results  have  led  a 
number  of  researchers  to  believe  that  the  set  of  30  test 
matrices  has  enough  variety  to  be  representati ve  of  typical 
application  problems  which  would  require  the  use  of  a  reduc¬ 
tion  algorithm.  An  interesting  example  throughout  testing 
reduction  algorithms  has  been  the  matrix  with  910  vertices, 
The  Snay  algorithm,  the  King  algorithm,  and  the  reverse  Cut- 
hill-Mcl<ee  algorithm  have  all  performed  poorly  on  it  indepen¬ 
dent  of  starting  points.  In  fact  there  is  no  starting  point 
for  which  Snay's  algorithm  or  King's  algorithm  performs  well 
on  this  example.  The  reason  why  is  an  open  problem. 
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